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Abstract. Motivated by the buckling of glassy crusts formed on evaporating 
droplets of polymer and colloid solutions, we numerically model the deformation 
and buckling of spherical elastic caps controlled by varying the volume between the 
shell and the substrate. This volume constraint mimics the incompressibility of 
the unevaporated solvent. Discontinuous buckling is found to occur for sufficiently 
thin and/or large contact angle shells, and robustly takes the form of a single 
circular region near the boundary that 'snaps' to an inverted shape, in contrast 
to externally pressurised shells. Scaling theory for shallow shells is shown to 
well approximate the critical buckling volume, the subsequent enlargement of the 
inverted region and the contact line force. 



PACS numbers: 46.32,+x, 46.70.De, 81.15.-z 

The properties of fluid interfaces are strongly modified by the presence of 
a permeable solid layer pp. This is deliberately exploited in solid-stabilised or 
'Pickering' emulsions, where added colliodal particles become pinned at the liquid- 
liquid interfaces, inhibiting dewetting and vastly reducing droplet coalescence [21 IS] ■ 
Solid interfacial layers may also form spontaneously, as in the observation of a glassy 
'crust' on the surface of evaporating polymer or colloid solutions in thin film 01 El, 
suspended droplet and partially-wetted droplet [7| 03 El El E3 geometries. In 
all of the above examples, it has been observed that, as the fluid volume(s) change 
with time, either by evaporation or draining (and might also be expected to occur 
during Ostwald ripening of emulsions), the interfacial area contracts, compressing the 
solid layer which then wrinkles or buckles as an elastic sheet. Predicting the onset and 
nature of this buckling is thus crucial to controlling system evolution and preventing 
the occurrence of undesired features in any given application. 

The deformation and buckling of thin elastic shells is a classic problem; see 
e.g. [131 1141 HTflllfij . However, the systems mentioned above introduce a complication 
that has not, to the best of our knowledge, been properly treated, namely that the 
presence of the incompressible fluid imposes a volume constraint on the space of 
allowed deformations (on time scales shorter than e.g. the evaporation time) [7]. It 
might be thought that controlling the volume V would be identical to imposing some 
corresponding pressure P, but this is not necesarily true with regards the stability 
of the shell: an S-shaped P-V curve would reveal different limit points, and hence 
distinct buckling events ^| ^j, depending on which quantity is being controlled 
(similar to S-shaped flow curves in non-linear rheology, which may be stable under 
stress control but unstable under an imposed flow rate or vice versa; see e.g. |19p. 
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Below we describe numerical simulations of a thin elastic shell with the same 
geometry as the crust on a partially- wetted spherical droplet, with the shell 
deformation driven by incrementally reducing the volume between it and the substrate. 
This geometry was chosen for its industrial relevance: the drying of inks and 
paints [201 IZ3 and 'bottom-up' manufacturing of microscale electronic devices via 
microfluidic technology [22] ah involve the evaporation of partially-wetted solutions. 
Our primary finding is of a single class of buckling event for thin and/or steep shells, in 
which a region of the shell near the fixed boundary 'snaps' to an inverted configuration 
[see Fig.^, in contrast to the range of axisymmetry preserving and breaking buckling 
modes already documented for the pressure-controlled case ^2]- For shallow shells, the 
critical buckling volume V c extracted from the simulations, as well as the subsequent 
enlargement of the inverted region, is shown to be well estimated by predictions of a 
simple scaling theory. We first describe the employed simulation method. 

Methodology: The crust is modelled as a thin elastic shell described by a 2-dimensional 
mid-surface S parameterised by a thickness h, such that the crust surfaces lie at a 
distance ±/i/2 normal to S. As the droplet surface is spherical, we take S to be a 
section of sphere of radius R that intersects the substrate at an angle 8o; see Fig.^c). 
Deviations from this unstressed state are described by a pair of rank- 2 tensors; Uij, the 
in-surface or membrane strains, and the change in curvature &Xij (with corrections 
for the difference between centroidal and mid-surfaces 221)- A standard closure due to 
Kirchoff and Love then gives the following expression for the strain energy, partitioned 
into a membrane (or 'stretch') contribution J7 s t re tchj and a flexural (or 'bend') term 



where E and v are respectively the Young's modulus and Poisson ratio of the shell 
material in bulk. Note that this is constitutivcly linear, all non-linearities are 
geometric in origin, arising from the finite displacements implicit in «y and A^y - 

Discretising JQ) with a regular mesh runs the risk of introducing undesired 
coupling between long-range mesh ordering and any buckling modes that may arise. 
To remove this possibility, we employ a statistically isotropic and homogeneous random 
mesh, formed by constructing a set of nodal points added with uniform probability 
on the spherical surface, with the constraint of a minimum distance between any two 
nodes. This was then Delaunay triangulated using qhull [23J - Energy contributions for 
local deformations are calculated in a simple finite element scheme, similar to that of 
recent full sphere |23] and wrinkling |2S] studies. Mesh edges are treated as Hookean 
springs of natural length Zo, so that 5C/ Strctch = \ot{l/h — I) 2 with I the current length, 
a = EhSA/(l — v 2 ) and 5 A the area assigned to this edge. The curvature between 
two adjacent faces, with unit normals hi and n 2 and whose midpoints are separated 
by a distance e [see Fig.^d)], is \ ~ 0/ e with (f> = arccos(fii ■ n 2 ) <C 1. The energy 
contribution is then SUbcnd = hP(x ~ Xo) 2 > where xo is the spontaneous curvature 
and (3 = h 2 a/12 to ensure convergence to the continuum limit suggested by QJ. 



E/bcnd d 
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Figure 1. Examples of (a) prebuckled and (b) post— buckled shell configurations 
(<?0 = 0-8, h/R = 0.03). (c) Undeformed geometry, a spherical cap of thickness h, 
contact angle 8o and radius of curvature _R. (d) Schematic of the discrete energy 
calculation (see text for details). 



A series of quasi-static configurations for a given mesh is generated by minimising 
the total strain energy for incrementally smaller volumes V by the non-linear 
conjugate gradient method, with the single scalar constraint V |23|- Clamped 
boundary conditions were employed, in which the position and angle of boundary 
elements are fixed, as suggested by experimental observation of real crusts 
Interactions with the substrate were introduced by including an energy cost U su b 
when shell nodes penetrate the substrate, with U su b made sufficiently stiff that the 
shell is essentially fiat when in contact, as confirmed by direct visualisation. 

Results: As the volume V is reduced from its undeformed value Vo, the shell 
initially undergoes axisymmetric deformations, with higher-order azimuthal modes 
(i.e. dimple — > mexican hat — > crater — > ...) arising for thinner and/or steeper shells, 
much as in the externally-pressurised case |28| . However, buckling events, identified 
here as discontinuous drops in the total strain energy for an arbitrarily small volume 
change, are qualitatively different, taking the form of a single, circular region near the 
boundary that becomes inverted, as seen in Fig. ^ (the proximity to the boundary is 
likely because it masks off part of the high-energy rim; see below). This has been 
observed even for contact angles 9q up to l\ radians, but differs from full spheres 
where multiple inverted regions are expected |26| and experimentally observed 

To quantify these observations, it is convenient to consider shallow shells 
SoCl, for which the two dimensionless parameters that define the geometry of the 
unperturbed shell, 9q and h/R, can be combined into the single scalar 



For externally pressurised shells, A is known to determine if a shell buckles, and the 
critical pressure at which it does ^3 EE]; we show here that it plays a similar role under 
volume control. Buckling only occurs for A > A c with Ac = 5.25 ±0.15, corresponding 
to thin and/or steep shells (note that this is significantly higher than the pressure 




(2) 



Volume-controlled buckling of spherical elastic caps 



4 



controlled value ,AP ressure « 3.3 |151 128| ). Furthermore, the dimensionless change in 
volume at which buckling occurs, AV c /Vb, approximately scales as A -2 as in Fig. [3 

This relation can be understood by deriving approximate scaling relations between 
various characteristic quantities, in the style of |29) . Let the characteristic radially- 
inward displacement be ( and assume this is the dominant mode. The undeformed 
height of the shell is #(l-cos0 o ) [see Fig.QJc)], so AV/Vb ~ (/R(l -cos0 o ) ~ (/R6 2 , 
for 6q <^1. Employing the known result that buckling occurs when £ ~ h [211123, we 
derive AV c /Vq ~ hjR9\ ~ A -2 , confirming the numerical findings. 

Extending this analysis to the strain energy provides insight into the physical 
meaning of the parameter A. Given the deformation mode described above, shell 
elements of area 5 A will be compressed by a characteristic membrane strain ~ 
C/R- For bending, note that, in the pre-buckled state, the shell deforms smoothly 
over the entire surface, i.e. over a length oc R9q- Since the curvature is related to 
the second derivative of transverse displacements, Ax^ ar ~ (/(R9q) 2 . Inspection 
of (JIJ then suggests characteristic energies SU stlc tch ~ Eh((/R) 2 SA and <5E/bcnd ~ 
Eh 3 [(/(R9 ) 2 ] 2 6A. The ratio is is SU stl - ct ch/SU hcnd ~ R 2 9%/h 2 - A 4 , so A > 1 
corresponds to a stretching dominated regime or membrane state, and conversely small 
A corresponds to a bending dominated response, somewhat comparable to flat plates. 

6 | , , 1 , 1 1 




X, = [12(1-v)] 1/4 (R9 2 /h) 1/2 



Figure 2. Critical volume scaled by A versus A in the buckling regime 
A > A c . rs 5.25 for the three different contact angles 8q shown in the legend. 
The horizontal line corresponds to AK ~ A -2 as predicted by scaling theory. 



In the post-buckled state, the inverted region has roughly the same radius of 
curvature as the undeformed state but opposite sign; this gives almost zero stretch 
energy, and is thus energetically favourable for small h. However, irrespective of h, 
the rim between inverted and non-inverted regions is highly curved and the bending 
stresses can never be ignored here. With this insight, it is possible to extend the scaling 
analysis to the post-buckled state. It has been shown that the characteristic width d of 
the rim that minimises the total energy scales as d ~ y/hR 29 . From purely geometric 
calculations, the volume of the inverted region of radius r is AV/Vb ~ ( r /R9o) 4 , thus 
scaling theory predicts d ~ \/h~R and r ~ R9 (AV/Vo) 1 ^ 4 in the post-buckled regime. 
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This is confirmed by the numerical data in Fig. EI Note that a parallel argument shows 
that the energy of n identical inverted regions scales as n 1 / 4 thus the robustly 
observed n — 1 case is indeed the energetically favoured one. 




Figure 3. Evolution of the buckled region as a function of the relative reduction 
in volume (Vb — V)/Vq = AV/Vb, found by fitting the strain energy density to the 
form Ce - " x— x ol _r ) / 2d after stereographic projection to the horizontal x-plane 
passing through the apex. The upper panel shows the radius r of the inverted 
region, along with with the 1/4 slope predicted by scaling theory. The lower panel 
gives the width d of the high-stress 'rim' scaled by the prediction d ~ V hR (the 
arrows indicate when the shell comes into contact with the base). 



In applications it is often important to know the forces at the contact line, as 
this can cause the crust to peel from the substrate, an undesirable effect in e.g. 
inks or paints. In the pre-buckled regime, the smoothly-varying membrane stresses 
scale with Ev%j ~ E^/R ~ EQqAV/Vq, so the force per length of contact line is 
/ ~ Eufj^h ~ Eh6^AV/Vo, which acts outwards since the shell is in compression. 
This is confirmed by the numerics; see Fig.H At the buckling volume AV C /V ~ h/R6l 
(when f c ~ Eh 2 /R), the contact force suddenly drops, as seen in the figure, where it 
remains of much lower magnitude until the shell becomes strongly deformed. Since 
this is now due to subdominant terms in the stress field, it is beyond the simple scaling 
theory presented here to predict its form. 

For completeness, we have also simulated a sample of hinged caps where the 
contact angle is free to change. Although the buckling is of a similar form, the critical 
volume obeys AV c /Vo ~ A~ Q with an exponent a as 2.25 rather than 2 as for the 
clamped case; also, the post-buckled configurations interact with the substrate much 
sooner, nullifying the simple post-buckling analysis presented above. 

Discussion: The broad conclusion of this work is the inequivalence of pressure and 
volume control with regards the stability of thin elastic shells. For the spherical 
cap geometry considered here, volume control reveals just a single class of instability, 
contrasting with externally-pressurised caps |15) . It is likely that a similar discrepancy 




Figure 4. The radial component of the contact line force per unit length 
f r = /cos(#o) ~ f versus Eh9^ AV/Vo for linear deformations, confirming the 
prediction of the scaling theory. (Inset) Non-linear behaviour for 6q = 0.4 and h 
decreasing in factors of 2 from h = 0.025 (top) to h = 0.0015625 (bottom). For 
all but the largest h, buckling occurs at the discontinuity in the curve. 



will arise in other geometries, suggesting established buckling criteria will need to be 
rethought whenever incompressibility introduces a volume constraint to the space of 
allowed deformation modes. 

Experiments on partially-wetted droplets have yet to show the buckling described 
here. This suggests that other mechanisms need to be incorporated, such as dynamics 
(solution and shell hydrodynamics and viscoelasticity), continued shell thickening 
during evaporation, and temperature gradients, amongst others. Indeed it is because of 
this the broad range of possibilities that a simple volume constraint was included in the 
simulations, rather than coupling to a uniform pressure field (the normal surface forces 
in the simulations are higher near the boundary); there seems little point demanding 
uniform pressure when the actual form will be far from trivial. It is hoped that a 
fuller model will soon be available to better understand this important phenomenon. 
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